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ABSTRACT 

Results  of  an  experimental  and  analytical  research  investigation  on  nozzle/afterbody 
drag  are  presented.  Experimental  afterbody  (and  boattail)  drag  coefficients  and  pressure 
distributions  are  discussed  for  an  isolated,  strut-mounted  nozzle/afterbody  model  for  the 
Mach  number  range  from  0.6  to  1.5.  Some  data  are  also  given  for  free-stream  unit  Reynolds 
numbers  from  one  million  to  approximately  four  million  per  foot.  The  experimental  data 
were  obtained  for  the  basic  model  with  an  air-cooled  and  a  water-cooled  Ethylene®/air 
combustor  to  provide  hot-jet  duplication  as  well  as  cold-jet  simulation.  The  temperature 
of  the  nozzle  exhaust  gas  was  varied  from  530°R  (burner  off)  to  approximately  2500°R 
for  several  nozzle  pressure  ratios  from  jet  off  to  those  corresponding  to  a  moderately 
under-expanded  exhaust  plume.  The  initial  series  of  experiments  was  conducted  with  the 
air-cooled  combustor,  and  the  effect  of  jet  temperature  on  afterbody  drag  was  somewhat 
masked  by  the  effects  of  the  secondary  airflow  from  the  cooling  air.  The  general  trend, 
however,  shows  a  decreasing  afterbody  drag  with  increasing  exhaust  gas  temperature  and 
with  decreasing  secondary  airflow  at  a  fixed  nozzle  pressure  ratio.  The  results  from  the 
water-cooled  combustor  show  similar  trends  with  increasing  exhaust  jet  temperature.  The 
differences  between  the  cold-jet  and  hot-jet  results  are  significant,  and  adjusting  the  cold-jet 
pressure  ratio  to  correct  for  the  changes  in  the  jet  specific  heat  ratio  with  temperature 
will  account  for  most  of  the  differences  observed.  A  parallel  analytical  investigation  is 
described  which  attempts  to  predict  the  afterbody  drag  of  isolated  nozzle/afterbody 
configurations.  The  effort  to  date  has  been  moderately  successful  in  accounting  for  the 
exhaust  plume  displacement  effects  (both  those  which  occur  because  of  the  inviscid 
expansion  and  those  which  occur  because  of  plume  mixing)  at  free-stream  conditions  where 
subsonic  flow  exists  over  the  entire  afterbody.  The  analytical  procedures  used,  with  their 
merits  and  shortcomings,  are  discussed. 
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SECTION  I 
INTRODUCTION 


A  significant  amount  of  testing  is  involved  in  developing  a  new  aircraft  and  power 
plant.  Traditionally,  airframe  tests  have  been  carried  out  with  careful  simulation  of  the 
external  flow  over  the  airframe  but  with  a  rather  approximate  simulation  of  the  flow 
issuing  from  the  engine.  Similarly,  engine  tests  have  been  conducted  with  very  careful 
simulation  of  the  flows  internal  to  the  engine  but  with  little  or  no  consideration  of  the 
external  flow. 

In  order  to  properly  evaluate  afterbody  drag,  realistic  simulation  of  the  internal  as 
well  as  the  external  flow  is  required.  The  hot-jet  exhaust  simulation  is  generally  obtained 
by  using  a  cold  fluid  (e.g.,  unheated  air)  and  adjusting  the  nozzle  pressure  ratio  to  match 
the  hot-jet  initial  inclination  angle.  Such  a  procedure  then  attempts  to  account  for  the 
differences  in  the  plume  expansion  caused  by  changes  in  the  jet  specific  heat  ratios  with 
exhaust  gas  temperature.  The  plume  displacement  effects  which  are  simulated  in  the  above 
manner  are  inviscid  in  nature,  and  although  in  theory  the  inviscid  plume  boundaries  are 
simulated,  the  effective  jet  boundaries  will  differ  because  of  the  viscous  mixing  at  the 
interface  of  the  exhaust  flow  and  the  external  boattail  flow.  The  nature  of  the  viscous 
region  depends  largely  on  the  gradients  existing  between  the  internal  nozzle  flow  and 
the  external  flow  over  the  afterbody.  To  investigate  this  viscous-inviscid  interaction,  which 
is  particularly  strong  at  transonic  speeds,  it  is  necessary  to  duplicate  rather  than  simulate 
the  engine  exhaust  gases.  The  research  reported  herein  is  a  parallel  experimental  and 
analytical  effort  to  determine  the  effects  of  the  hot  jet  on  the  boattail  drag.  The  analytical 
investigation  is  a  continuing  effort  to  predict  the  afterbody  drag  on  nozzle/afterbody 
models.  These  predictions  will  have  to  ultimately  account  for  the  effect  of  jet  entrainment 
of  the  boattail  flow  as  well  as  the  plume  displacements  caused  by  viscous  mixing  on  the 
afterbody  drag.  With  sufficiently  accurate  prediction  techniques  available,  the  proper 
procedure  for  simulating  a  hot-jet  exhaust  flow  with  a  cold  fluid  may  be  developed. 


The  experimental  portion  of  the  research  reported  herein  was  conducted  in  two  phases, 
which  were  differentiated  by  the  combustor  design.  Phase  I  was  conducted  using  an 
air/ethylenc  combustor  which  had  air  cooling  around  the  combustor  liner.  During  the  Phase 

I  investigation,  it  was  observed  that  the  secondary  airflow  (cooling'air)  affected  the  boattail 
pressure  distribution,  hence  producing  drag.  This  effect  opposed  the  trends  resulting  from 
the  introduction  of  the  hot  jet.  Because  of  the  secondary  airflow  effect,  the  water-cooled 
combustor  was  designed  (Phase  II)  to  eliminate  secondary  airflow  and  isolate  the  jet  exhaust 
temperature  effect  on  boattail  and  afterbody  drag.  During  Phase  I,  two  nozzle 
configurations  representing  a  typical  turbojet  nozzle  installation  and  a  typical  turbofan 
with  low-bypass  (1:1)  nozzle  installation  at  a  military  power  setting  were  used.  The  Phase 

II  investigation  utilized  the  turbojet  nozzle  installation  at  a  military  power  setting.  Data 
were  obtained  for  these  configurations  over  the  transonic  Mach  number  range. 


1 


AEDC-TR-74-9 


SECTION  II 
APPARATUS 


2.1  WIND  TUNNEL 

The  AEDC  16-Ft  Propulsion  Wind  Tunnel  (16T)  is  a  continuous  flow,  closed-circuit 
wind  tunnel  capable  of  operation  within  a  Mach  number  range  from  0.20  to  1 .60.  Tunnel 
16T  can  be  operated  within  a  stagnation  pressure  range  from  120  to  4000  psfa,  depending 
on  Mach  number,  with  a  stagnation  temperature  variation  capability  from  approximately 
80°F  to  a  maximum  of  160°F.  Tunnel  air  is  removed  and  replaced  with  conditioned  make-up 
air  from  an  atmospheric  dryer  to  facilitate  control  of  vitiation  caused  by  combustion  and 
to  control  the  specific  humidity  of  the  tunnel  air. 

2.2  EXPERIMENTAL  HARDWARE 

The  experimental  hardware  used  during  the  investigation  reported  herein  was  the  Air 
Force  Flight  Dynamics  Laboratory  (AFFDL)/Lockheed  isolated  nozzle  model  previously 
used  for  cold-flow  nozzle  studies.  The  model,  described  fully  in  Ref.  1,  was  modified  to 
accommodate  an  Ethylene®/air  combustor,  which  was  used  to  provide  hot  exhaust  jet 
duplication  of  a  typical  turbojet  operation  (Section  2.2.1).  The  AFFDL/ Lockheed  isolated 
nozzle  model  is  an  axisymmetric  body  with  an  overall  length  of  approximately  1 53  in.  The 
model  had  a  14-deg,  half-angle  conical  nose,  which  was  faired  into  the  primary  model 
diameter  of  9.86  in.  A  boundary-layer  trip  consisting  of  0.055-in.-diam  steel  spheres 
^  spot-welded  to  a  trigging  at  a  circumferential  spacing  of  four  sphere  diameters  was  located 
on  the  conical  nose'  1 2  In.  aft  of  the  cone  vertex.  The  external  surfaces  of  the  afterbody  and 
boattail  used~_wereuthe  AFFDL/Lockheed  configuration  CDE1  (Convergent-Divergent 
Ejector,  military  power  setting)  described  in  Ref.  2.  The  model  was  mounted  in  Tunnel  16T 
on  a  strut  with  an  aft  sweep  angle  of  31.8  deg.  A  sketch  showing  the  basic  external  model 
dimensions  and  the  model  location  in  the  wind  tunnel  test  section  is  shown  in  Fig.  1 
(Appendix  I).  Installation  photographs  of  the  model  are  presented  in  Fig.  2. 

The  internal  model  configurations  deviated  from  the  AFFDL/Lockheed  CDE1 
configuration  in  order  to  accommodate  the  combustor.  High-pressure  air,  which  was  used 
to  simulate  the  nozzle  exhaust  for  the  cold-flow  portion  of  testing  and  which  was  used 
as.  an  oxidizer  in  the  combustion  of  Ethylene,  was  ducted  through  the  strut  to  a  plenum 
in  the  forward  section  of  the  model. 

2.2.1  Combustor 

An  Ethylene/air  combustor  based  on  a  National  Advisory  Committee  for  Aeronautics 
(NACA)  design  (Ref.  2)  was  used  to  provide  hot  exhaust  jet  duplication.  (Ethylene  (C2  H4 ) 
is  a  gaseous  hydrocarbon  fuel  which,  when  burned  in  air,  produces  exhaust  products  which 
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very  closely  duplicate  the  exhaust  products  of  JP-4  burned  in  air.)  Two  combustor 
configurations  were  used  during  the  experimental  investigation.  The  overall  dimensions 
of  the  two  configurations  were  similar,  and  many  components  were  interchangeable.  The 
first  combustor  used  secondary  airflow  for  cooling  the  combustor  liner  and  nozzle,  whereas 
the  second  combustor  and  nozzle  were  water  cooled. 

The  flow  straighteners  forward  of  the  combustor  and  the  flameholder  were  the  same 
for  both  combustors.  The  flameholder  consisted  of  four  doughnut  rings  interconnected 
with  fuel  flow  passages,  and  the  forward  doughnut  ring  served  as  the  fuel  injection  ring. 
Ignition  was  accomplished  by  injecting  a  small  quantity  of  tri-ethyl  borane  (TEB,  a 
pyrophoric  fuel)  into  a  retainer  mounted  on  the  flameholder.  Flow  mixers  were  placed 
at  the  rear  of  the  flameholder  at  radial  positions  of  45,  135,  225,  and  315  deg  (0  deg 
represented  the  top  of  the  model)  to  mix  the  hot  flow  from  the  combustion  zone  with 
the  air  around  the  outside  of  the  flameholders.  This  was  done  for  both  combustors  to 
provide  a  uniform  temperature  profile  at  the  exit  of  the  combustor.  The  physical  differences 
in  these  configurations  because  of  the  different  cooling  schemes  is  explained  in  the 
following  paragraphs.  A  schematic  of  the  basic  burner  and  a  photograph  of  the  flameholder 
are  shown  in  Fig.  3. 

2.2.1. 1  Air-Cooled  Combustor 

The  air-cooled  combustor  liner  was  constructed  from  stainless  steel  and  had  an  overall 
length  of  27.45  in.,  an  inside  diameter  of  i6.15  in.,  and  a  wall  thickness  of  0.125  in. 
(Fig.  4a).  The  combustor  liner  terminated  in  a  sonic  throat  with  a  contraction  ratio  of 
2.36.  Secondary  airflow  for  cooling  was  ducted  to  the  secondary  flow  annulus  from  the 
primary  stream  through  a  variable  position  sleeve  valve  (Fig.  4a).  The  two  air  streams, 
primary  and  secondary,  were  mixed  at  the  exit  plane  of  the  combustor  sonic  nozzle  and 
were  exited  through  the  common  nozzle.  Flow  conditioners  in  both  the  primary  and  the 
secondary  air  streams  were  used  to  provide  a  uniform  flow  field  to  the  combustor  and 
secondary  flow  annulus,  respectively. 

2.2.1. 2  Water-Cooled  Combustor 

The  water-cooled  combustor  was  constructed  from  copper  with  an  inside  diameter 
of  6.15  in.  and  a  wall  thickness  of  0.125  in.  (Fig.  4b).  The  combustor  terminated  with 
a  convergent-divergent  nozzle  with  a  throat  (diameter  of  3.6  in.  The  contraction  ratio  for 
this  combustor  was  2.92.  Cooling  water  was  supplied  to  the  combustor  through  the  model 
strut  and  flow  passages  internal  to  the  model  (Fig.  4b).  A  continuous  water  flow  rate  of  70 
gal/min  was  maintained  throughout  the  investigation. 

2.2.2  Nozzle  Configurations 

Two  nozzle  configurations  were  used  in  conjunction  with  the  air-cooled  combustor. 
Each  configuration  utilized  the  AFFDL/ Lockheed  convergent-divergent  ejector  afterbody 
and  boattail  external  surfaces  (Configuration  CDE1)  and  force  balances.  The  internal 
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military  nozzle  configuration  was  constructed  from  stainless  steel  and  had  a  throat  diameter 
of  4.8  in.  and  an  expansion  ratio  of  1.0,  with  a  length-to-diameter  ratio  of  1.25  (Fig. 
5a).  Since  the  burner  was  operated  choked,  the  effective  aTea  ratio  of  the  nozzle  was 
1.44  for  this  configuration.  A  secondary  air  passage  sleeve  was  used  to  reduce  the  annulus 
area  around  the  combustor,  thereby  restricting  the  secondary  airflow  to  just  the  amount 
required  for  cooling  the  combustor  liner  (approximately  15  percent  of  primary  airflow). 
The  secondary  air  passage  sleeve  was  removed  for  the  1:1  bypass  nozzle  configuration, 
allowing  a  secondary  airflow  rate  nearly  equal  to  the  primary  stream.  The  1:1  bypass 
nozzle  configuration  was  also  constructed  of  stainless  steel  and  had  an  expansion  ratio 
of  1.06  (Fig.  5b).  The  nozzle  was  preceded  by  a  7.5-in.  mixing  chamber  to  better  simulate 
the  mixing  region  of  the  hot  core  and  cold  annulus  air  experienced  in  turbofan  engines. 

The  internal  nozzle  configuration  used  with  the  water-cooled  combustor  was  a 
convergent-divergent  (CD)  military  nozzle  with  a  nominal  design  area  ratio  of  1.43  and  a 
divergence  angle  of  6.4  deg.  Because  of  design  requirements  imposed  by  water  cooling,  the 
base  area  (4.78  in.2)  of  the  CD  nozzle  was  significantly  larger  than  the  air-cooled  military 
turbojet  nozzle  (Fig.  5c). 

2.3  INSTRUMENTATION 

Critical  flow  venturi  meters  were  used  to  measure  the  total  air  and  Ethylene  flows. 
Secondary  airflow  rate  for  the  air-cooled  combustor  was  determined  from  pressure  and 
temperature  measurements  in  the  secondary  flow  annulus  at  a  known  cross-sectional  area. 
The  primary  airflow  was  then  found  by  subtracting  the  secondary  airflow  from  the  total 
airflow  through  the  venturi. 

Two  thin-walled  cylindrical  ring  force  balances  were  used  on  the  air-cooled  combustor 
configuration  for  measurement  of  nozzle  external  drag  and  nacelle  afterbody  drag, 
respectively.  The  balances  contained  metallic  bellows  seals  that  prevented  air  from  leaking 
through  the  balance  flexures.  For  each  of  the  two  cylindrical  ring  balances,  the  effective 
area  of  the  bellows,  which  was  determined  from  balance  calibration,  was  used  to  make 
appropriate  pressure-area  corrections  to  the  balance  data.  Static  pressures  along  the  body 
external  surfaces  were  also  measured. 

The  combustion  chamber  pressure  probe  was  located  adjacent  to  the  tri-ethyl  borane 
line  (Fig.  3).  The  military  nozzle  air-cooled  combustor  configuration  had  six  static  pressure 
orifices  at  approximately  the  15-  and  200-deg  radial  locations  (two  rows  of  three  orifices 
located  internally,  approximately  180  deg  apart).  The  1:1  bypass  nozzle  had  only  four 
nozzle  static  pressure  orifices,  with  an  additional  total  and  static  pressure  measurement 
made  in  the  secondary  flow  passage.  The  location  of  the  internal  pressure  instrumentation 
for  the  two  air-cooled  combustor  nozzles  is  shown  in  Figs.  5a  and  b.  In  addition,  the 
water-cooled  combustor  configuration  was  instrumented  with  three  static  pressure  orifices 
in  the  combustor  and  two  static  pressure  orifices  in  the  divergent  nozzle  (Fig.  5c). 

Measurement  uncertainty  values  are  presented  in  Table  I  (Appendix  II)  for  some  of 
the-  more  important  parameters. 
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SECTION  III 
PROCEDURES 

3.1  EXPERIMENTAL  PROCEDURES 

The  experimental  investigation  was  divided  into  two  phases.  Phase  I  utilized  the 
air-cooled  combustor  configuration  and  the  two  associated  nozzles,  whereas  Phase  II  utilized 
the  water-cooled  combustor  configuration.  The  experimental  procedures  were  slightly 
different  for  these  two  phases  and  will  be  enumerated  separately. 

Phase  I  Procedure 

1.  Wind  tunnel  Mach  number  and  Reynolds  number  conditions  were 
established,  and  jet-off  data  were  obtained. 

2.  Cold  flow  was  established  through  the  model  at  a  specified  setting  of  the 
secondary  airflow  rotating  sleeve  valve,  and  a  series  of  nozzle  pressure  ratios 
were  set.  Data  were  obtained  at  each  ratio. 

3.  The  total  airflow  through  the  model  was  fixed,  and  the  rotating  sleeve  valve 
was  varied  from  its  fully  closed  to  its  fully  open  position  to  obtain  secondary 
airflow  effects.  This  was  usually  done  at  a  nozzle  pressure  ratio  near  that 
set  for  hot-flow  operation. 

4.  The  combustor  mixture  was  ignited  and,  for  a  given  fuel/air  ratio,  a  nozzle 
static  pressure  ratio  was  set  (using  pe3/p„  as  the  controlling  ratio).  Data 
were  obtained  at  four  basic  fuel/air  ratios  as  near  to  the  same  nozzle  pressure 
ratio  as  possible. 

When  setting  the  selected  pressure  ratio  with  the  hot  flow,  both  the  fuel  and  the  airflow 
were  varied  at  a  given  fuel/air  ratio. 

Phase  II  Procedure 

1.  Wind  tunnel  Mach  number  and  Reynolds  number  conditions  were 
established,  and  jet-off  data  were  obtained. 

2.  Cold  flow  was  established  through  the  model,  and  a  series  of  nozzle  pressure 
ratios  were  set.  Data  were  obtained  at  each  ratio. 

3.  The  combustor  mixture  was  ignited,  and  the  fuel/air  ratio  was  set  to  obtain 
the  desired  temperature.  A  series  of  nozzle  pressure  ratios  were  set  by 
varying  fuel  and  airflow  to  maintain  a  constant  fuel/air  ratio.  Data  were 
obtained '  at  each  nozzle  pressure  ratio.  The  nozzle  pressure  ratio  survey 
was  conducted  at  four  fuel/air  ratios. 
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Schlieren  photographs  were  obtained  at  each  of  the  nozzle  pressure  ratios  for  the 
hot  exhaust  jet  duplication  and  at  selected  pressure  ratios  for  the  cold  exhaust  jet 
simulation. 

3.2  DATA  REDUCTION  PROCEDURES 


The  primary  performance  parameters  presented  in  this  report  are  the  afterbody  and 
boattail  drag  coefficients  obtained  from  an  integration  of  the  experimental  pressure 
distributions  and  the  computed  skin  friction.  Both  drag  coefficients  were  calculated  using 
similar  relationships.  The  relationships  used  in  calculating  the  boattail  drag  coefficient  are 
presented  here  to  demonstrate  the  general  technique  employed. 


Cdbt  -  Cdbtp  +  Cdbtf 
The  pressure  drag  coefficient  is 

144 (DBTpp) 

CdBTP - qls 


(1) 


(2) 


where  q„  is  the  free-stream  dynamic  pressure  and  S  is  the  projected  area  based  on  the 
maximum  diameter  of  the  model.  The  boattail  pressure  drag  (DBTpp)  is  calculated  from 
the  experimental  pressures  by  the  following  equation: 

i=448 

Dbtpp  =  0.21 048 6p„  -  ^  PiSi  (3) 

i=4  03 

Here,  i  represents  an  individual  pressure  measurement  location'  (see  Tables  II  and  III); 
Pi  is  the  experimental  pressure,  and  Sj  is  the  projected  area  dver  which  the  measured 
pressure  acts.  The  friction  drag  coefficient  is  calculated  from  the  relationship 

144(DBTf) 

Cdbtf  =  - ; -  (4) 

q«S 

where  Dbtf  is  the  total  friction  force  obtained  by  summation  of  the  incremental  friction 
forces 

i=4  4  8 

Dbtf  =  ^  (Cfb)  (q«)  (Aj)  (cos  0j)  (5) 

i=  4  0  3 

where  Cfb  is  the  average  skin-friction  coefficient  based  on  the  local  Reynolds  number, 
qg  is  the  local  dynamic  pressure,  Aj  is  the  area  of  the  external  surface  segment,  and  04 
is  the  local  boattail  surface  angle,  The  afterbody  friction  drag  coefficient  is  obtained  from 
the  empirical  correlation  of  local  skin-friction  coefficient  versus  Reynolds  number 
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developed  by  Sivells  and  Payne  (Ref.  3)  for  incompressible  turbulent  flow  over  a  flat 
plate  with  zero  pressure  gradient. 

0.088  (log  Rei  -  2.3686) 

Cfi  =  (log  Rei  -  1.5)3 

The  internal  performance  parameters  of  the  exhaust  jet  were  obtained  in  the  following 
manner  for  the  water-cooled  and  air-cooled  models. 

In  the  case  of  the  water-cooled  model,  three  static  pressure  orifices  located  at  various 
axial  positions  in  the  water-cooled  combustion  chamber  (Fig.  5c)  were  used  to  obtain 
the  chamber  static  pressure.  The  pressure  measured  by  the  orifice  farthest  downstream 
Gust  ahead  of  the  start  of  the  throat  contraction)  was  used  to  calculate  the  exhaust  jet 
total  pressure,  ptj,  in  the  following  manner: 


Ptj  -  Pc*  ^  p  ^  (7) 

where  pt/p  is  the  subsonic  isentropic  pressure  ratio  based  on  the  chamber  contraction 
ratio  of  2.92.  The  isentropic  pressure  ratio  is  a  weak  function  of  the  ratio  of  specific 
heats  (pt/p  decreases  2.0  percent  for  a  change  in  7  from  1.40  to  1.30).  Because  of  this 
dependence  on  7,  ptj  was  calculated  with  an  iteration  procedure  wherein  the  initial 
calculation  was  made  using  an  assumed  exhaust  jet  temperature.  Exhaust  jet  temperature 
was  calculated  in  the  iteration  loop  using  the  relationship  given  in  Eq.  (8).  The  physical 
geometry  of  the  chamber  and  the  measured  pressure  and  mass  flow  rate  were  used  in 
the  calculation 


Tg 


as 


Ptj  Kf  At  k.z 

wt 


(8) 


where  Kp  is  Fleigner's  constant  and  is  a  function  of  7,  Aj  is  the  physical  throat  area, 
Kz  is  the  nozzle  throat  discharge  coefficient  determined  from  cold-flow  test  of  the  nozzle, 
and  wt  is  the  total  mass  flow  rate  in  the  chamber.  The  total  mass  flow  is  given  by 

wt  =  wa  +  Wf  (9) 

where  wa  is  the  chamber  airflow  rate  measured  by  a  critical  flow  venturi  and  Wf  is 
the  ethylene  flow  rate  calculated  from  the  Tg  a  s  value  in  the  iteration  loop  and  an  assumed 
combustion  efficiency  of  80  percent.  The  entire  iteration  procedure  was  made  necessary 
because  a  change  in  the  thermodynamic  quality  of  the  gaseous  ethylene  (from  a  perfect 
gas  to  a  saturated  vapor)  resulted  in  a  deviation  from  the  perfect  gas  law  and  prevented 
the  venturi  measuring  system  (originally  designed  to  measure  a  gaseous  ethylene  flow) 
from  giving  accurate  flow  measurements. 
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In  the  case  of  the  air-cooled  model  tests,  the  problem  encountered  with  ethylene 
quality  was  also  suspected,  which  made  the  gas  temperature  an  unknown  since  fuel  flow 
was  not  known.  In  this  model  the  only  pressure  measured  in  the  combustion  chamber 
was  made  with  the  extra,  open  TEB  supply  line  shown  in  Fig.  3a.  This  orifice  did  not 
give  an  accurate  chamber  pressure  with  cold  flow  in  the  chamber  because  of  a  relatively 
high  air  velocity  in  that  region  of  the  chamber.  However,  when  there  was  combustion 
in  the  chamber  a  correlation  of  the  ptJ-  calculated  for  the  water-cooled  model  with  the/ 
pressure  measured  with  the  TEB  line  showed  less  than  a  1 -percent  deviation  between' 
the  two.  Based  on  this  correlation,  Eq.  (8)  was  again  used  to  compute  the  actual  gas 
temperature  for  the  hot  flow  conditions  using  the  pressure  from  the  TEB  line  for  ptj. 
The  main  uncertainty  in  this  procedure  was  the  discharge  coefficient  for  the  air-cooled 
nozzle,  which  was  assumed  equal  to  unity. 

The  complexity  of  the  combustor  and  nozzle  fabrication  resulted  in  the  loss  of  some  of 
the  nozzle  pressure  taps  which  were  planned.  Because  of  the  loss  of  these  pressure  taps,  two 
additional  pressures  were  calculated  for  use  in  presenting  the  data:  a  calculated  nozzle  exit 
static  pressure,  pe,  for  the  water-cooled  model  and  a  measured  average  nozzle  exit  static 
pressure,  peav,  for  the  air-cooled  model.  The  water-cooled  pe  was  calculated  from  Eq.  (10) 
as  follows: 

P.  =  P.J  (~)  00) 

where  p/ptJ-  is  the  iseritropic  pressure  ratio  based  on  the  effective  area  ratio  of  the  nozzle 
(1.48)  and  is  a  function  of  y.  This  calculated  nozzle  exit  pressure  was  higher  than  the 
measured  pressure  (Pes)  at  the  nozzle  exit  by  approximately  2  percent.  [However,  the 
measured  p£5  was  probably  reduced  slightly  by  the  low  pressure  in  the  base  region  of 
the  nozzle.  The  nozzle  exit  static  pressure  for  the  air-cooled  combustor  configuration 
(Peav)  was  taken  to  be  the  arithmetic  average  (Eq.  11)  of  the  two  measured  pressures 
180  deg  apart  located  near  the  nozzle  exit  plane  (Fig.  5a).  Thus, 

1 

Peav  =  J  (PE 3  +  PE6)  (11) 

SECTION  IV 

EXPERIMENTAL  RESULTS 


4.1  GENERAL 

The  experimental  investigation  was  conducted  in  two  phases.  Phase  I  utilized  the 
air-cooled  combustor  configuration  with  both  a  military  turbojet  and  a  1:1  bypass  military 
turbofan  nozzle  configuration.  Data  were  obtained  using  the  turbojet  nozzle  configuration 
at  Mach  numbers  of  0.6,  0.9,  1.1,  1.2,  and  1.5  over  the  nozzle  pressure  ratio  range  shown 
in  Fig.  6  at  a  Reynolds  number  of  2.5  x  1 06 /ft.  Four  nominal  temperatures  (1000®, 
1500°,  2000®,  and  2500"R)  were  investigated  at  each  Mach  number.  Data  were  also 
obtained  with  the  turbojet  nozzle  using  cold  exhaust  jet  simulation  for  comparison  with 
the  hot  exhaust  jet  duplication. 
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The  1 : 1  bypass  turbofan  nozzle  configuation  was  also  investigated  at  Mach  numbers 
of  0.6,  0.9,  1.1,  1.2,  and  1.5  at  a  Reynolds  number  of  2.5  x  1 06 /ft.  Two  primary  stream 
temperatures  of  approximately  1000°  and  1500°R  were  investigated.  Cold  exhaust  jet 
simulation  data  for  comparison  with  the  hot-jet  duplication  conditions  were  also  obtained. 
Cold-  and  hot-plume  data  (at  the  specified  fuel/air  ratios)  were  also  obtained  at  a  Reynolds 
number  of  4.75  x  106/ft  and  Mach  numbers  of  0.6,  0.9,  and  1.1. 

Phase  II  of  the  experimental  investigation  was  conducted  using  the  water-cooled 
combustor  configuration.  The  investigation  was  conducted  at  a  Reynolds  number  of  2.5 
x  1 06  /ft  over  the  same  Mach  number  range  as  Phase  I,  at  four  nominal  exhaust  gas 
temperatures.  The  cold  exhaust  jet  simulation  conditions  were  repeated  during  Phase  II 
for  data  comparison  with  the  hot  exhaust  jet  duplication.  During  Phase  II,  an  extensive 
Reynolds  number  survey  (from  1.0  x  1 06  /ft  to  3.0  x  1 06  /ft)  was  also  conducted  at  Mach 
numbers  of  0.6  and  0.9  over  the  nominat-fuel/air  ratio  range.  Schlieren  photographs  were 
taken  at  selective  nozzle  pressure  ratios  to  supplement  the  pressure  data  obtained-  during 
Phase  II. 

The  afterbody  and  boattail  drag  coefficients  were  chosen  as  the  primary  parameters 
to  compare  nozzle/ afterbody  performance  with  hot  exhaust'  jet  duplication  to  the 
performance  obtained  using  cold  exhaust  plume  simulation.  Two  methods  of  determining 
these  drag  coefficients  were  used  during  Phase  I.  The  first  method  utilized  the  force 
balances  and  measured  drag  force,  whereas  the  second  method  obtained  drag  by  an 
integration  of  measured  pressure  distribution  over  the  afterbody  and  boattail  with  a 
computed  skin  friction  added.  A  comparison  of  the  boattail  drag  coefficients  obtained 
using  each  method  is  shown  in  Fig.  7.  The  agreement  between  the  two  methods  is  good. 
Because  of  the  close  agreement  with  balance  data  and  the  ever-present  possibility  of  balance 
shifts  caused  by  the  internal  complexity  of  the  model,  the  drag  coefficients  obtained 
by  pressure  integration  plus  a  computed  skin  friction  are  presented  in  this  report. 

4.2  WATER-COOLED  COMBUSTOR  CONFIGURATION  (PHASE  II) 

Typical  pressure  distributions  over  the  boattail  obtained  from  the  water-cooled 
combustor  are  shown  in  Figs.  8a  through  d.  The  effect  that  changing  the  exhaust  jet 
nozzle  pressure  ratio  has  on  the  boattail  pressure  distribution  is  shown  in  Figs.  8a  and 
b  for  Mach  numbers  0.9  and  1.1,  respectively.  For  both  Mach  numbers,  increasing  the 
nozzle  pressure  ratio  caused  an  increase  in  the  pressure  on  the  boattail.  The  effect  was 
felt  upstream  of  the  boattail  on  the  afterbody  at  Mach  number  0.9  as  shown  by  the 
pressure  difference  at  the  forwardmost  station  of  the  boattail.  At  Mach  1.1,  an  imbedded 
shock  occurred  on  the  boattail  and  restricted  the  effect  of  the  exhaust  plume  to  the 
portion  of  the  boattail  downstream  of  the  shock.  However,  the  shock  tended  to  move 
upstream  with  increasing  nozzle  pressure  ratio.  The  effect  of  exhaust  plume  temperature 
on  pressure  distribution  at  Mach  numbers  0.9  and  1.1  is  shown  in  Figs.  8c  and  d.  Increasing 
exhaust  temperature  generally  causes  an  increase  in  pressure  on  the  boattail  at  Mach 
number  0.9,  with  the  effect  being  felt  on  the  afterbody  as  well.  At  Mach  1 .1 ,  the  imbedded 
shock  limited  the  effect  of  temperature  to  the  segment  of  the  boattail  downstream  of 


9 


AEDC -TR-74-9 


the  shock  location.  To  illustrate  the  extent  of  the  influence  of  the  exhaust  plume,  the 
pressure  distribution  along  essentially  the  entire  model  length  is  shown  in  Fig.  8e  for 
the  conditions  given  in  Fig.  8a.  Schlieren  photographs  at  the  three  exhaust  plume 
temperatures  shown  in  Fig.  8d  at  Mach  number  1.1  are  presented  in  Figs.  9a  through 
c. 


The  component  drag  coefficients  (boattail  drag  coefficient,  CdbTj  afterbody  drag 
coefficient,  Cd  ab  .  and  total  drag  coefficient,  Cdt -which  is  the  sum  of  Cdbt  and  Cdab) 
as  a  function  of  nozzle  pressure  ratio  for  each  of  the  five  Mach  numbers  and  a  Reynolds 
number  of  2.5  x  106/ft  are  shown  in  Figs.  10a  through  e.  In  each  case,  the  cold-plume 
simulation  data  are  presented  with  the  varying  hot-plume  duplication  for  comparison. 
Above  a  nozzle  pressure  ratio  of  5.66  (nozzle  completely  filled),  Cdbt  decreased  with 
increasing  pressure  ratio.  The  level  of  drag  coefficient  obtained  with  the  hot-exhaust 
duplication  at  a  given  nozzle  pressure  ratio  generally  decreased  with  increasing  exhaust 
temperature.  At  the  subsonic  Mach  numbers  0.6  and  0.9,  the  afterbody  drag  coefficient, 
Cd  ab  .  decreased  slightly  with  both  pressure  ratio  and  exhaust  plume  temperature  increases 
(Figs.  10a  and  b).  At  the  supersonic  Mach  numbers  of  1.1,  1.2,  and  1.5,  however,  Cdab 
was  completely  insensitive  to  changes  in  either  of  the  variables.  This  insensitivity  to  pressure 
ratio  and  exhaust  plume  temperature  changes  is  caused  by  the  inability  of  the  pressure 
ratio  and  temperature  effects  to  move  upstream  of  the  imbedded  shock  on  the  boattail 
(Figs.  8b  and  d). 

The  mechanism  which  governed  the  behavior  of  boattail  drag  was  investigated  by 
Bergman  in  Ref.  4.  The  qualitative  results  of  his  investigation  are  presented  in  Fig.  1 1 . 
Bergman  proposed  that  the  boattail  drag  coefficient'  (cold  flow)  was  dominated  by  two 
physical  effects.  The  first  was  a  displacement  outward  of  the  external  flow  streamlines 
caused  by  the  exhaust  plume  emitting  from  the  model  nozzle.  The  outward  displacement 
of  the  streamlines  resulted  in  stronger  flow  recomprqssion 1  on  the  boattail  and  thereby 
had  a  beneficial  effect  on  drag.  The  shape  of  the  exhaust  plume  is  a  function  of  nozzle 
pressure  ratio  for  fixed  tunnel  conditions.  The  second  effect,  entrainment,  began  when 
the  exhaust  plume  velocity  was  approximately  equal  to  free-stream  velocity.  It  increased 
boattail  drag  by  increasing  boattail  velocity  and  effectively  lowering  the  boattail  pressure. 
The  net  effect  on  boattail  drag  is  a  superposition  of  the  two  effects. 

The  cold-flow  Cdbt  obtained  during  this  investigation  had  the  general  shape  and 
characteristics  obtained  by  Bergman  (Fig.  10a).  No  data  were  obtained  at  the  low  nozzle 
pressure  ratios  (less  than  4.0)  with  a  hot  exhaust;  however,  it  is  hypothesized  that  the 
general  behavior  will  be  retained  at  the  higher  temperatures.  For  a  fixed  pressure  ratio, 
the  hot  exhaust  plume  will  spread  more  than  the  cold  plume,  causing  a  lower  drag  on 
the  boattail.  The  effect  of  exhaust  plume  temperature  on  entrainment  is  unknown, 
however,  and  the  entire  problem  is  the  subject  of  a  continuing  analytical  investigation 
guided  by  the  experimental  data  obtained  during  this  investigation. 
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4.2.1  Matching  of  the  Plume  Inclination  Angle 

If  the  changes  in  boattail  drag  caused  by  changes  in  nozzle  pressure  ratio  are  primarily 
inviscid  in  nature  (i.e.,  caused  by  changes  in  the  inviscid  plume  shape),  then  corrections 
to  the  cold  flow  nozzle  pressure  ratio  can  be  made  to  simulate  the  change  in  plume 
shape  produced  by  changes  in  the  exhaust  gas  temperature.  The  effect  of  increasing  the 
exhaust  gas  temperature  is  to  decrease  the  specific  heat  ratio  (7)  of  the  exhaust  gas  with 
a  resulting  increase  in  the  exhaust  plume  initial  inclination  angle  and  maximum  diameter. 
Thus,  in  order  to  simulate  a  hot-jet  plume  (7  <  1.4)  operating  at  a  given  pressure  ratio 
using  cold  air  (7  =  1 .4),  the  nozzle  pressure  ratio  must  be  increased  to  offset  the  effect 
of  7.  The  correction  for  7  effects  is  to  adjust  the  nozzle  pressure  ratio  at  a  constant 
value  of  drag  coefficient  to  match  initial  jet  inclination  angles.  The  measured  performance 
with  the  theoretical  change  in  the  nozzle  pressure  ratio  predicted  from  changes  in  7  at 
the  five  free-stream  Mach  numbers  investigated  are  compared  in  Figs.  12a  through  e.  The 
correction  for  7  compensates  for  a  large  percentage  of  the  exhaust  gas  temperature-induced 
effects.  The  effect  is  corrected  more  accurately  at  the  higher  nozzle  pressure  ratios  than 
at  the  lower  pressure  ratios.  The  effectiveness  of  the  7  correction  also  appears  to  be 
dependent  on  the  free-stream  Mach  number.  The  7  correction  agrees  more  closely  at 
M„  =  1.2.  At  free-stream  Mach  numbers  less  than  1.2  the  effect  of  plume  temperature 
is  slightly  larger  than  the  7  effects  predict,  while  at  M„  =  1.5,  the  7  correction  is  larger 
than  the  measured  effects  of  plume  temperature. 

4.2.2  Method  of  Presenting  Data 

Historically,  nozzle/afterbody  drag  has  been  presented  as  a  function  of  total  nozzle 
pressure  ratio  (NPR),  defined  as  Ptj/p„.  For  a  fixed  area  ratio  nozzle  this  introduces  into 
consideration  a  varying  nozzle  exit  static  to  total  pressure  ratio  (pe/ptj)  with  7  which 
must  be  properly  accounted  for  when  the  hot-jet  effects  on  afterbody  and  boattail 
performance  are  evaluated.  This  variation  of  pe/ptj  with  7  can  be  eliminated  from  the 
observed  effects  if  the  nozzle/afterbody  performance  is  presented  as  a  function  of  nozzle 
static  pressure  ratio  (pc/p«).  A  comparison  of  the  two  methods  of  data  presentation  is 
shown  in  Figs.  13a  and  b  for  total  drag  coefficient,  Cox.  at  Mach  1.1.  At  a  constant 
value  of  NPR  the  effect  of  jet  temperature  on  drag  appears  to  be  more  pronounced  than 
it  is  for  constant  values  of  pc/p«.  This  is  caused  by  the  variation  in  pe/Pt j  with  7  discussed 
above.  The  manner  of  variation  of  pe /Pt j  *s  shown  by  superimposing  on  Fig.  13a  lines 
of  constant  pc/p».  A  presentation  of  the  exhaust  jet  temperature  effects  on 
nozzle/afterbody  performance  as  a  function  of  pe/p»  (Fig.  13b)  tends  to  collapse  the 
data  spread  observed  in  Fig.  13a. 

The  method  of  data  presentation  in  this  report  is  dictated  by  the  measured  parameters 
available.  In  the  Phase  II  investigation  the  total  nozzle  pressure  ratio,  NPR,  was  chosen 
because  the  measured  chamber  pressure,  ptj,  was  more  reliable  than  was  the  nozzle  exit 
static  pressure  measurement.  For  the  Phase  1  investigation,  the  nozzle  exit  static  pressure 
measurement  was  more  reliable,  and  therefore  the  data  were  presented  as  a  function  of 
nozzle  exit  static  pressure  ratio.  Both  methods  are  valid  and  present  the  same  information 
as  long  as  the  variation  of  pe/p»  with  7  is  considered. 
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4.2.3  Reynolds  Number  Effect 

To  determine  the  effect  of  Reynolds  number  on  boattail  drag  coefficient,  a  Reynolds 
number  survey  from  1.0  x  106  to  3.0  x  1 06  /ft  was  conducted  at  Mach  numbers  of  0.6 
and  0.9.  Drag  data  were  obtained  with  cold  and  hot  exhaust  plumes.  The  effect  of  Reynolds 
number  on  boattail  pressure  distribution  is  presented  in  Fig.  14  at  Mach  0.6  for  the 
cold-flow  data.  The  effect  of  varying  Reynolds  number  from  1.5  x  106  to  3.0  x  106 /ft 
is  imperceptible.  The  variation  of  the  boattail  pressure  drag  with  Reynolds  number  at 
a  constant  nozzle  pressure  ratio  and  various  gas  temperatures  is  shown  in  Fig.  15.  At 
both  Mach  0.6  and  0.9,  the  effect  of  Reynolds  number  on  Cdt  is  very  small  over  the 
range  investigated.  Figures  16a  through  f  present  the  drag  coefficients  obtained  at  Mach 
0.6  and  0.9  at  each  of  the  Reynolds  numbers  as  a  function  of  nozzle  pressure  ratio. 

4.3  AIR-COOLED  COMBUSTOR  CONFIGURATION  (PHASE  I) 

Phases  I  and  II  of  the  experimental  investigation  were  conducted  using  the  same 
external  afterbody  and  boattail  configuration  (Section  2.2.2).  The  size  of  the  base  region 
was  enlarged  for  the  Phase  II  model  because  of  requirements  imposed  by  water  cooling 
the  combustor  and  nozzle  (Fig.  5).  This  larger  base  area  affected  the  level  of  the  boattail 
drag,  resulting  in  greater  drag  at  a  given  nozzle  static  pressure  ratio  than  was  measured 
on  the  Phase  I  model  (Fig.  17).  A  correction  for  base  area  effect  for  the  cold  exhaust 
plume  at  Mach  0.9  was  made  using  the  results  obtained  in  Ref.  5  and  shown  in  Fig. 
17  to  illustrate  that  the  base  area  was  responsible  for  the  Phase  II  increased  drag.  The 
corrected  Phase  II  drag  coefficient  was  comparable  to  the  level  measured  in  Phase  I. 
Sufficient  base  pressure  data  to  allow  correction  of  Cdbt  f°r  the  remainder  of  the  data 
were  not  available;  therefore,  in  reviewing  the  results  of  Phase  I  and  comparing  them 
to  Phase  II,  one  should  remember  that  the  effect  of  base  area  was  present  and  the  absolute 
levels  will  not  be  comparable.  The  trends  remain  the  same,  however,  for  both  phases. 

4.3.1  Effect  of  Secondary  Airflow  on  Nozzle  Performance  (Phase  I) 

Secondary  airflow  to  provide  combustor  cooling  during  hot  operation  for  the 
air-cooled  combustor  was  introducted  into  the  nozzle  as  described  in  Section  2.2.  During 
cold-flow  operation,  the  secondary  airflow  was  maintained,  and  it  affected  the  internal 
nozzle  performance  as  shown  in  Figs.  18  and  19.  The  nozzle  pressure  ratio  as  a  function 
of  axial  station  increased  with  increasing  Secondary  airflow  ratio  (W2/WO,  indicating  that 
a  larger  annular  portion  of  the  available  nozzle  area  was  used  by  the  lower-energy, 
secondary  airflow.  During  an  analysis  of  the  experimental  data,  a  disagreement  between 
the  two  nozzle  static  pressures  (pe3  and  pe6)  was  observed.  These  pressures  are  located 
in  a  common  axial  plane  1  in.  upstream  of  the  nozzle  exit  and  approximately  180  deg 
apart.  Throughout  the  cold-exhaust  plume  investigation,  the  agreement  between  these 
pressures  was  affected  by  W2/W1  (Fig.  19).  Agreement  was  good  for  airflow  ratios  between 
four  and  ten  percent;  however,  for  airflow  ratios  less  than  four  or  greater  than  ten  percent, 
Pe3  and  pe6  differed  significantly.  This  disagreement  has  not  been  satisfactorily  explained, 
but  it  was  probably  a  combination  of  secondary  airflow  effects  and  combustor 
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misalignment.  This  internal  nozzle  effect  was  transmitted  to  the  boattail  external  forces 
through  the  change  in  plume  shape  and  size.  Figures  20a  through  e  are  a  presentation 
of  the  boattail  drag  coefficient  (cold  flow)  as  a  function  of  measured  nozzle  static  pressure 
ratio.  Lines  of  constant  W2/W1  are  identified.  As  W2/W1  increased,  boattail  drag  coefficient 
increased  at  a  constant  nozzle  static  pressure  ratio  of  all  Mach  numbers,  with  the  exception 
of  Mach  0.6.  This  effect  of  secondary  airflow  on  the  boattail  performance  partially  obscured 
the  effects  of  exhaust-jet  temperature. 

4.3.2  Effect  of  Hot  Exhaust  on  Nozzle  Performance  (Phase  I) 

A  presentation  of  the  component  drag  coefficients  as  a  function  of  nozzle  static 
pressure  ratio  for  each  of  five  Mach  numbers  is  shown  in  Figs.  21a  through  e.  The  drag 
coefficients  were  obtained  over  a  wide  range  of  exhaust  gas  temperatures,  using  the 
cold-exhaust  plume  to  establish  the  relationship  between  boattail  drag  coefficient  and 
pressure  ratio.  Because  of  the  disagreement  between  the  measured  values  of  pe3  and  pgg, 
an  average  of  the  two  pressures  was  used  to  represent  the  nozzle  static  pressure  ratio. 
The  boattail  drag  coefficient  for  the  hot-exhaust  plume  was  then  determined  at  specific 
values  of  nozzle  static  pressure  ratio  and  was  compared  to  the  cold-plume  drag  coefficient 
on  the  assumption  that  the  relationship  between  drag  coefficient  and  pressure  ratio  remains 
generally  the  same  even  though  the  exhaust  gas  temperature  is  increased.  The  experimental 
results  of  Phase  II  verified  this  assumption.  The  effects  of  exhaust  gas  temperature  and 
increasing  secondary  airflow  have  opposing  trends  and  tend  to  cloud  the  effect  of  exhaust 
plume  temperature. 

4.3.3  One-To-One  (1:1)  Bypass  Nozzle 

Component  drag  coefficients  measured  on  the  1 : 1  bypass  nozzle  configuration  at 
each  of  the  Mach  numbers  (0.6,  0.9,  1.1,  1.2,  and  1.5)  are  presented  in  Figs.  22a  through 
e  as  a  function  of  nozzle  pressure  ratio.  The  component  drag  coefficients  were  determined 
at  a  unit  Reynolds  number  of  4.75  x  106  for  the  1:1  bypass  nozzle  configuration  at 
Mach  numbers  of  0.6,  0.9,  and  1.1  and  are  compared  in  Figs.  22a  through  c  with  the 
coefficients  obtained  at  2.5  x  1 06  /ft.  Both  Cdbt  and  Cdab  at  Mach  numbers  0.6  and 
0.9  were  higher  at  a  Reynolds  number  of  4.75  x  10®  /ft  than  at  2.5  x  1 06  /ft.  The  effect 
of  Reynolds  number  was  more  pronounced  on  Cdab  than  on  Cdbt-  At  a  free-stream 
Mach  number  of  1.1,  the  Reynolds  number  change  had  no  effect  on  Cdbt  and  only 
a  slight  effect  on  Cdab- 


SECTION  V 

THEORETICAL  ANALYSIS 

5.1  PHYSICAL  STATEMENT  OF  PROBLEM 

The  initial  attempts  to  provide  analytical  prediction  of  the  pressure  distribution  and 
drag  of  an  isolated  nozzle/afterbody  configuration  were  carried  out  under  the  usual 
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assumption  that  the  viscous-inviscid  interaction  could  be  accounted  for  by  an  iteration 
scheme.  Such  a  scheme  involves  an  iteration  between  the  five  distinct  but  dependent  regions 
of  the  total  flow  field  which  are  identified  in  Fig.  23.  Region  A  represents  the  inviscid 
core  of  the  internal  nozzle  and  is  normally  supersonic.  Region  B,  also  assumed  to  be 
inviscid  in  nature,  is  the  external  flow  around  the  model  and  can  be  subsonic,  supersonic, 
or  mixed  flow.  Regions  C  and  D  are  viscous  boundary  layers  internal  to  the  nozzle  and 
external  to  the  body,  respectively.  At  the  lip  of  the  nozzle,  the  viscous  boundary  layers 
merge  to  form  a  viscous  mixing  region,  E.  In  general,  the  boundary  conditions  for  each 
region  are  provided  by  the  solution  of  the  adjacent  region.  The  inviscid  regions,  A  and 
B,  can  be  defined  by  specifying  either  a  solid  boundary  or  a  pressure  boundary.  Region 
A  will  be  defined  by  a  solid  boundary  consisting  of  the  nozzle  plus  the  displacement 
thickness  from  Region  C  and  by  the  external  displacement  boundary  from  Region  E. 
Region  B  will  be  defined  by  a  solid  boundary  consisting  of  the  boattail  plus  the 
displacement  thickness  from  Region  D  and  by  the  internal  displacement  boundary  from 
Region  E.  The  solutions  of  Regions  A  and  B  will  provide  the  pressures  at  the  boundaries. 
These  pressures,  plus  the  physical  nozzle  and  afterbody  coordinates,  allow  a  solution  of 
the  boundary  layers,  Regions  C  and  D.  Conditions  from  Regions  C  and  D  are  used  as 
starting  conditions  for  the  viscous  mixing  region,  E.  The  axial  pressure  gradient  and  the 
location  of  a  reference  line  are  also  required  for  Region  E.  The  internal  and  external 
displacement  thicknesses  are  computed  with  respect  to  this  reference  line.  The  boundary 
conditions  for  each  region  listed  are  summarized  in  Table  IV. 

5.2  NUMERICAL  METHODS 

The  computational  technique  used  to  solve  the  descriptive  equations  for  the  entire 
flow  field  involves  an  iteration  between  the  inviscid  flow  fields,  A  and  B,  and  the  viscous 
flow  fields,  C,  D,  and  E.  The  procedure  currently  used  for  carrying  out  the  iteration 
is  given  below: 

1.  The  nozzle  flow  properties  in  Region  A  are  computed  using  the  nozzle 
physical  coordinates  and  a  quiescent  atmosphere  at  the  free-stream  static 
pressure  in  Region  B. 

2.  The  physical  coordinates  of  the  afterbody  and  the  inviscid  boundary  of 
Region  A  are  used  to  compute  the  flow-field  properties  in  Region  B. 

3.  The  inviscid  flow  field  properties  in  Regions  A  and  B  are  matched  by 
recomputing  the  plume  expansion  (Region  A)  with  the  computed  pressure 
gradient  from  Region  B.  A  single  iteration  is  usually  sufficient  for  this 
matching. 

4.  The  turbulent,  viscous  boundary-layer  characteristics  in  Regions  C  and  D 
are  computed  using  the  pressure  distributions  from  Regions  A  and  B, 
respectively,  along  with  the  physical  coordinates  of  the  model. 
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5.  The  flow  properties  in  the  free-shear  turbulent  mixing  layer  in  Region  E 
are  computed  using  the  inputs  from  the  viscous  boundary  layers  C  and 
D  and  the  inviscid  reference  line  from  Region  A. 

6.  The  displacements  and  boundary-layer  characteristics  from  the  viscous  layers 
(Regions  C,  D,  and  E)  are  used  to  correct  the  model  plus  the  plume 
boundary  to  recompute  the  flow  fields  A  and  B  to  begin  the  second 
iteration.  Region  A  may  or  may  not  need  to  be  recomputed  since  small 
corrections  are  usually  produced  on  this  flow  from  the  viscous  effects. 

Computer  programs  have  been  developed  for  the  numerical  solution  of  the  flow 
properties  in  the  regions  discussed  (Table  V).  Region  A  can  be  solved  using  the 
axisymmetric  Method  of  Characteristics  Program  developed  by  Lockheed  (Ref.  6).  This 
program  is  used  for  supersonic  flow  and  is  also  employed  to  solve  Region  B  when  the 
external  flow  field  is  supersonic.  If  Region  B  is  subsonic,  a  program  developed  by  Pratt 
&  Whitney  (Ref.  7)  is  used  for  its  solution.  This  solution  is  carried  out  using  a  numerical 
solution  of  the  nonlinear,  transonic  flow  equation.  The  equation  is  written  in  the  form 
of  a  Poisson  equation  with  the  nonlinear  terms  on  the  left-hand  side.  The  solution  is 
then  performed  by  the  method  of  successive  substitutions  starting  with  the  solution  of 
the  Laplace  equation.  The  treatment  of  imbedded  shocks  limits  the  practical  utility  of 
this  program  to  afterbody  flows  which  are  subsonic  everywhere.  The  turbulent 
boundary-layer  calcuations,  over  either  the  external  afterbody/boattail  or  the  internal 
nozzle  (Regions  C  and  D),  are  made  with  the  implicit  finite-difference  scheme  of  Patankar 
(Ref.  8).  The  viscous  mixing  region  E  can  also  be  solved  with  the  Patankar  program, 
employing  a  second  option  using  boundary-layer  assumptions. 

The  scope  of  this  report  does  not  allow  a  detailed  description  of  the  calculation 
techniques  used  in  these  various  computer  programs.  However,  the  reader  is  directed  to 
the  cited  references  for  a  more  complete  description. 

The  computational  technique  outlined  above  does  not  include  the  effect  of 
entrainment  on  the  boattail  pressure  distribution.  This  effect  is  necessary  for  a  complete 
analytical  solution  to  the  problem,  and  efforts  to  incorporate  the  entrainment  effect  on 
the  boattail  pressures  in  the  computational  technique  are  being  carried  out  at  AEDC.  The 
computational  technique  is  also  inadequate  when  nearly  normal  shockwaves  are  located 
on  the  boattail  of  the  afterbody. 

5.3  COMPARISON  OF  THEORY  AND  EXPERIMENTAL  RESULTS 

The  methods  described  were  used  to  compute  the  pressure  fields  over  the  aft  portion 
of  the  isolated  nozzle  model  and,  finally,  the  drag  over  the  afterbody  and  boattail.  These 
efforts  are  reasonably  successful  at  the  subsonic  Mach  numbers  for  the  air-cooled  model, 
which  has  a  small  base  area.  During  the  course  of  malting  the  drag  calculations,  it  was 
found  that  certain  approximations  could  be  made  to  shorten  the  iteration  process  when 
the  external  stream  was  subsonic.  First,  because  of  the  low  nozzle  pressure  ratios  used 
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during  subsonic  flight  (see  Fig.  6),  matching  the  nozzle  exit  static  pressure  to  free-stream 
static  pressure  ratio  is  the  primary  parameter  required  to  match  the  inviscid  plume  shapes. 
Hence,  changing  the  plume  total  temperature  from  520  to  2700°R  had  only  a  negligible 
effect  on  the  inviscid  plume  boundaries,  which  are  shown  in  Figs.  24  and  25  for  M„ 
=  0.6  and  0.9,  respectively.  Secondly,  mixing  was  primarily  from  the  high-energy  plume 
into  the  low-energy  external  flow,  and  the  internal  displacement  boundary  did  not  vary 
significantly  with  respect  to  temperature.  These  two  effects  combine  to  allow  the  external 
flow  and  boundary-layer-mixing  iteration  to  be  carried  out  without  affecting  the  basic 
plume  boundary  calculation. 

Although  the  inviscid  pressure  distribution  could  be  calculated  initially  without  any 
boundary-layer  displacement  corrections  to  the  afterbody/boattail  contour,  an  attempt  was 
made  to  shorten  the  iteration  cycle  by  making  the  displacement  corrections  from 
boundary-layer  calculations  carried  out  using  the  experimental  pressure  distribution.  The 
resulting  displacement  thickness  for  M*,  =  0.6  is  shown  as  the  dashed  line  in  Fig.  24. 
The  external  displacement  boundaries  for  four  plume  temperatures  are  also  shown.  Initial 
conditions  for  the  viscous  mixing  region  were  matched  to  the  internal  and  external 
boundary-layer  characteristics  at  the  nozzle  exit  plane.  Using  these  results  as  an  input 
boundary  for  the  inviscid  external  flow  solution,  a  much  steeper  pressure  gradient  was 
computed  on  the  boattail  than  was  measured.  When  this  predicted  pressure  was  imposed 
on  the  external  boundary  layer,  separation  was  predicted  at  the  point  noted  in  Fig.  24. 

As  can  be  seen  from  this  result,  the  attempt  to  shorten  the  viscous-inviscid  iteration 
cycle  by  using  the  boundary-layer  characteristics  from  the  experimental  pressure  was  not 
successful  since  the  second  iteration  indicated  boundary-layer  separation  not  predicted  from 
the  experimental  pressure  distribution.  The  calculation  procedure  outlined  in  Section  5.2 
does  not  have  any  method  available  to  predict  displacement  thicknesses  in  the  separated 
region.  This  requires  a  detailed  prediction  of  the  reverse  flow  profiles  and  the  extent  of 
the  separation  (or  abnormally  thick  boundary  layer)  region.  Recently  there  has  come  to 
the  attention  of  the  authors  a  calculation  procedure  developed  for  the  separated  turbulent 
boundary  layer  imbedded  in  the  transonic  flow  over  a  bump  on  a  wind  tunnel  wall  (Ref. 
9).  It  is  possible  that  such  a  procedure  can  be  incorporated  into  the  method  outlined 
previously,  and  an  attempt  will  be  made  to  do  so.  Since  such  a  calculation  procedure 
was  not  available  at  the  time  of  writing,  it  was  necessary  to  smoothly  fair  in  an  effective 
displacement  boundary  from  the  point  of  separation  to  the  effective  plume  displacement 
boundary  in  order  to  predict  the  experimental  pressure  distribution  from  the  corrected 
body  shape.  The  fairing  was  arbitrary  and  was  carried  out  in  a  manner  so  that  the 
experimental  data  agreed  with  the  analytical  results.  Efforts  are  still  under  way  to  develop 
reasonable  criteria  to  use  in  making  the  correction  to  the  body  and  plume  displacement. 
Two  points  appear  evident  when  the  displacement  corrections  are  made:  (1)  the 
downstream  variation  of  the  initial  plume  displacement  effects  remain  very  small  and  can 
be  predicted  during  the  first  cycle  of  the  viscous-inviscid  iteration,  and  (2)  the  changes 
in  the  pressure  distribution  between  iterations  cause  very  little  change  in  the  boundary-layer 
characteristics  ahead  of  the  point  where  separation  was  initially  predicted.  This  small  change 
was  probably  caused  by  (1)  the  weak  viscous-inviscid  iteration  present  on  the  particular 
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afterbody  used  and  (2)  the  fact  that  the  analysis  was  done  at  subsonic  Mach  numbers 
where  shock/boundary-layer  iterations  are  not  present.  jThe  strong  iteraction  region  from 
the  apparent  separation  point  to  the  apparent  point  of  attachment  to  the  exhaust  plume 
does  have  to  be  iterated  on  to  predict  the  correct  experimental  afterbody  pressure 
distribution.  The  final  fairings  for  the  displacement  boundaries  used  in  the  external  flow 
calculations  are  shown  in  Fig.  24  for  M„,  =  0.6.  Similar  results  are  shown  in  Fig.  25  for 
M»  =  0.9.  Preliminary  calculations  indicated  that  a  region  of  locally  supersonic  flow  existed 
over  the  shoulder  of  the  afterbody  at  M_  =  0.9.  The  afterbody  shape  was  modified 
somewhat  to  eliminate  this  locally  supersonic  region  so  that  the  external  flow  could  be 
treated  as  wholly  subsonic.  It  was  realized  that  the  pressures  would  not  be  well  predicted 
in  the  modified  region;  however,  it  was  felt  that  these  effects  of  modification  would  be 
sufficiently  localized  so  that  the  pressure  near  the  end  of  the  boattail  would  be  adequately 
predicted. 

Using  the  boundaries  shown  in  Figs.  24  and  25,  external  flow  solutions  were  carried 
out  for  =  0.6  and  0.9.  Pressure  distributions  over  the  body  are  shown  in  Fig.  26 
at  =  0.6  for  different  nozzle  plume  temperatures.  The  numerical  method  predicts  the 
general  data  trend  over  the  body.  It  is  noted  that  the  method  predicts  boattail  pressures 
that  increase  at  an  increasing  rate  with  temperature.  The  data  generally  support  this  trend, 
although  the  agreement  is  not  perfect. 

Pressure  distributions  over  the  body  are  shown  in  Fig.  27  for  =  0.9.  Similar 
theoretical  trends  with  temperature  and  agreement  with  the  data  over  the  boattail  were 
noted  for  M,*,  =  0.6.  The  pressure  over  the  shoulder  of  the  afterbody  is  not  well  predicted 
in  the  region  from  10  to  16  in.  upstream  of  the  nozzle  exit,  since  the  body  was  modified 
in  this  area  for  the  analytical  model. 

The  pressure  and  skin  friction  over  the  boattail  shown  in  Figs.  26  and  27  were 
integrated  to  obtain  drag.  The  results  of  these  integrations  are  presented  in  Fig.  28  as 
boattail  drag  coefficient  versus  exhaust  gas  temperature  and  are  compared  with  the 
experimental  data  obtained  with  both  the  air-cooled  and  the  water-cooled  combustors. 
The  slope  of  the  curve  of  drag  coefficient  as  a  function  of  exhaust  gas  temperature  for 
the  two  experimental  configurations  agrees  closely,  while  the  slope  obtained  from  theory 
I  is  steeper.  The  absolute  levels  of  drag  coefficient  predicted  by  theory  and  obtained 
experimentally  with  the  air-cooled  combustor  configuration  are  in  fair  agreement  at  an 
exhaust  gas  temperature  of  530’R.  The  absolute  level  of  the  drag  coefficient  obtained 
with  the  water-cooled  combustor  at  530°R  exhaust  gas  temperature  was  higher  than  those 
obtained  from  the  air-cooled  combustor  and  theoretical  predictions.  This  difference  is 
attributed  to  the  large  base  area  of  the  water-cooled  combustor  discussed  previously. 

A  plot  of  the  pressure  data  shown  in  Fig.  29  at  M»  =  1.2  is  typical  of  the  data 
when  the  external  stream  is  supersonic.  A  Method  of  Characteristics  calculation  is  shown 
for  comparison  in  Fig.  29.  The  data  and  the  theoretical  prediction  agree  quite  well  over 
the  shoulder  of  the  afterbody,  but  not  over  the  boattail.  First,  there  was  apparently  a 
disturbance  at  the  metric  break  between  the  afterbody  and  the  boattail  (6  in.  upstream 
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of  the  nozzle  exit)  that  was  not  predicted  by  the  MOC.  Second,  a  shock  was  located 
approximately  4  in.  upstream  of  the  nozzle  exit.  There  is  no  method  currently  available 
to  predict  either  the  location  or  the  strength  of  this  shock.  To  assume  that  the  flow 
separates  immediately  behind  the  shock  and  reattaches  near  the  maximum  diameter  of 
the  plume  would  require  a  turning  angle  which  yields  a  pressure  rise  across  the  shock 
of  the  same  order  as  was  measured.  However,  this  type  of  analysis  is  not  sufficiently 
accurate  to  evaluate  the  drag  over  the  afterbody. 

SECTION  VI 
CONCLUSIONS 

The  effect  of  exhaust  plume  temperature  on  afterbody  and  boattail  drag  was 
investigated  for  two  combustor  and  three  internal  nozzle  configurations  installed  in  an 
isolated  pod.  The  investigation  provided  both  experimental  and  analytical  data  for  the 
internal  and  external  flow  fields.  The  significant  results  and  conclusions  are  summarized 
below.  (It  should  be  emphasized  that  the  following  results  and  conclusions  were  obtained 
for  a  particular  isolated  external  boattail  configuration  and  may  not  necessarily  be  the 
same  for  other  configurations.) 

1.  For  the  water-cooled  combustor  configuration  with  a  turbojet  military 
nozzle,  the  maximum  change  in  total  drag  coefficient  at  a  Reynolds 
number  of  2.5  x  1 06 /ft  was  160  boattail  drag  counts  at  Mach  1.1  and 
90  counts  at  Mach  0.6  over  an  exhaust  gas  temperature  range  from  530° 
to  approximately  2700°R.  For  each  of  the  Mach  numbers  investigated,  drag 
coefficient  decreased  monotonically  with  increasing  exhaust  gas  temperature 
at  a  constant  nozzle  pressure  ratio. 

2.  From  the  hot-flow  data  obtained,  it  appears  that  correcting  the  cold-flow 
exhaust  pressure  ratio  for  changes  in  specific  heat  ratio  will  compensate  for  a 
large  precentage  of  the  effect  of  a  hot  exhaust  on  afterbody  performance. 

The  apparent  displacement  and  entrainment  effects  produced  by  differences 
in  viscous  mixing  between  hot  and  cold  exhausts  also  have  to  be  considered. 

3.  Exhaust  plume  temperature  affected  the  afterbody  drag  coefficient  at  Mach 
0.6  and  0.9.  The  effect  was  negligible  at  Mach  1.1,  1.2,  and  1.5. 

4.  A  Reynolds  number  survey  from  1.0  x  106  to  3.0  x  1 06 /ft  was  conducted 
at  Mach  0.6  and  0.9  for  the  water-cooled  combustor  configuration.  The 
effect  of  Reynolds  number  variation  on  total  boattail  drag  coefficient  was 
less  than  30  drag  counts  over  the  range  of  the  investigation  for  both  Mach 
numbers  0.6  and  0.9. 
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5.  The  boattail  drag  coefficient  obtained  with  the  air-cooled  combustor, 
military  nozzle  configuration  during  the  cold-exhaust  plume  investigation 
showed  a  strong  influence  of  secondary  airflow.  When  presented  as  a 
function  of  average  nozzle  static  pressure  ratio,  the  boattail  drag  increased 
with  an  increase  in  percentage  of  secondary  airflow. 

6.  The  boattail  drag  coefficient  obtained  with  the  air-cooled  combustor/ 
military  nozzle  configuration  behaved  in  the  same  manner  as  that 
obtained  with  the  water-cooled  configuration.  The  maximum  change  in 
boattail  drag  coefficient  at  a  Reynolds  number  of  2.5  x  106/ft  was 
182  drag  counts  at  Mach  0.9  and  50  drag  counts  at  Mach  1.5  over 
an  exhaust  gas  temperature  range  from  530°  to  approximately  2900°R. 

7.  The  boattail  drag  coefficient  decreased  with  increasing  exhaust  jet 
temperature  at  all  of  the  free-stream  Mach  numbers  of  the  investigation 
for  the  1:1  bypass  turbofan  nozzle  configuration.  The  temperature  effect 
agreed  in  general  with  the  results  obtained  with  the  water-cooled/ military 
nozzle  configuration. 

8.  Analytical  techniques  to  calculate  the  effect  of  temperature  on  boattail  drag 
coefficient  were  begun  with  emphasis  placed  on  subsonic  Mach  numbers 
0.6  and  0.9.  The  developed  method  agrees  well  with  the  experimental  results 
obtained  from  the  air-cooled  combustor.  This  method  does  not  account, 
however,  for  the  effect  of  entrainment  and  base  area  on  boattail  pressure 
distribution. 
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a.  View  Looking  Downstream 
Fig.  2  Model  Installation  Photographs 
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b.  View  Looking  Upstream 
Fig.  2  Concluded 
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a.  Sketch  of  Combustor 
Fig.  3  Combustor  Details 
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b.  Flameholder 
Fig.  3  Concluded 
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AFFDL/LOCKHEED  CONVERGENT- DIVERGENT  EJECTOR  NOZZLE 


a.  Air-Cooled  Combustor 
Fig.  4  Model  Internal  Details 
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AFFDUlockheed  Convergent-Divergent  Ejector 


b.  Water-Cooled  Combustor 
Fig.  4  Concluded 
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a.  Air-Cooled  Combustor,  Turbojet  Military  Nozzle 
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b.  Air-Cooled  Combustor,  Turbofan  Military  Nozzle  (1:1  Bypass  Ratio) 


c.  Water-Cooled  Combustor,  Turbojet  Military  Nozzle 
Fig.  5  Schematic  of  Nozzle  Configurations  and  Nozzle  Pressure  Instrumentation 
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Fig.  7  Comparison  of  Pressure  Integration  Boattail  Drag 
Coefficient  with  Force  Balance  Boattail  Drag 
Coefficient  at  M„,  =  1.2  and  Re  =  2.5  x  1 06  /ft 
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AXIAL  DISTANCE  FROM  NOSE.  in. 

a.  Effect  of  Static  Pressure  Ratio  at  -  0.9 
Fig.  8  Typical  Boattail  Static  Pressure  Distributions 
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e.  Pressure  Distribution  over  Model  at  M„  =  0.9 
Fig.  8  Concluded 
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a.  TG as  =  530° R,  NPR  =  9.17 

Fig.  9  Typical  Schlieren  Photographs  at  =  1.1  and  Re  =  2.5  x  1 06 /ft 
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b.  Tq As  =  1182°R,  NPR  =  9.07 
Fig.  9  Continued 
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C.  Tgas  =  2464° R,  NPR  =  9.10 
Fig.  9  Concluded 
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a.  =  0.6 

Fig.  10  Effect  of  Nozzle  Pressure  Ratio  on  Turbojet  Nozzle  Component 
Drag  Coefficients  at  Various  Exhaust  Jet  Temperatures  Using  a 
Water-Cooled  Combustor  at  Re  =  2.5  x  106/ft 
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Fig.  11  Qualitative  Effect  of  Entrainment  and  Plume  Shape  on  Boattail  Drag 
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Fig.  12  Concluded 
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a.  Boattail  Drag  Coefficient  as  a  Function  of 
Nozzle  Total  Pressure  Ratio  at  M_  =  1.1 


0.4  0.6  1.2  1.6  2.0  2.4 

Pe/p- 


b.  Boattail  Drag  Coefficient  as  a  Function  of 
Nozzle  Static  Pressure  Ratio  at  M„  =  1 .1 
Fig.  13  Comparison  of  Two  Methods  of  Presenting  Drag  Coefficient  Data 
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Fig.  14  Effect  of  Reynolds  Number  on  Boattail  Pressure  Distribution  at  -  0.6 
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a.  M„  =  0.6 


b.  NL  =  0.9 

Fig.  15  Effect  of  Reynolds  Number  on  Boattail  Pressure  Drag  Coefficient  at  Various 
Exhaust  Jet  Temperatures  Using  a  Water-Cooled  Combustor  at  a  Nozzle 
Pressure  Ratio  of  1.5 
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b.  IVL  =  0.6,  Re  =  2.0  x  106/ft 

Fig.  16  Effect  of  Nozzle  Pressure  Ratio  on  Component  Drag  Coefficients  at 
Various  Reynolds  Numbers  Using  a  Water-Cooled  Combustor 
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d.  M„  =  0.6,  Re  =  3.0  x  106/ft 
Fig.  16  Continued 
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M.  =  0.9,  Re  =  1.0  x  106/ft 
Fig.  16  Continued 
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Fig.  17  Comparison  of  Boattail  Drag  Coefficient  Levels  Using  Water-Cooled 
and  Air-Cooled  Combustors  at  M„  =  0.9 
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Fig.  18  Military  Nozzle  Internal  Pressure  Distribution  Cold  Flow  Using  an 
Air-Cooled  Combustor 
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Fig.  21  Effect  of  Nozzle  Pressure  Ratio  on  Turbojet  Nozzle  Component  Drag 
Coefficients  at  Various  Exhaust  Jet  Temperatures  Using  an  Air-Cooled 
Combustor  at  Re  =  2.5  x  10B/ft 
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b.  IVL  =  0.9 
Fig.  21  Continued 
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Fig.  21  Continued 
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Fig.  21  Continued 
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a.  fVL  =  0.6 

■  Fig.  22  Effect  of  Nozzle  Pressure  Ratio  on  Turbofan  Nozzle  Component  Drag 
Coefficients  at  Various  Exhaust  Jet  Temperatures  Using  an  Air-Cooled 
Combustor 
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Fig.  22  Continued 
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Fig.  22  Continued 


IVL=  1.2,  Re  =  2.5  x  10* /ft 
Fig.  22  Continued 


Fig.  23  Schematic  of  the  Various  Flow  Regions  involved  in  the  Afterbody  Drag  Calculation  Problem 
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Fig.  24  Effective  Displacement  Boundaries  for  the  External  Flow  at  M.  =  0.6  and  Re  =  2.5  x  IQS/ft 
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Fig.  25  Effective  Displacement  Boundaries  for  the  External  Flow  at  =  0.9  and  Re  =  2.5  x  10®/ft 
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Fig.  26  Afterbody  Pressure  Distribution:  Comparison  Between  Experiment  and  Theory  for  Various 
Plume  Temperatures  at  M„  =  0.6  and  Re  =  2.5  x  106 /ft 
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Fig.  27  Afterbody  Pressure  Distribution:  Comparison  Between  Experiment  and  Theory  for  Various 
Plume  Temperatures  at  M«,  =  0.9  and  Re  =  2.5  x  106/ft 
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O  AIR  COOLED  COMBUSTOR 
A  WATER  COOLED  COMBUSTOR 
□  THEORY 


EXHAUST  GAS  TEMPERATURE, T6AS , °R 

a.  M„  =  0.6 


EXHAUST  GAS  TEMPERATURE,  T6AS ,  *R 


b.  =  0.9 

Fig.  28.  Comparison  of  Theoretical  and  Experimental  Boattail  Drag 
Coefficients  at  Re  =  2.5  x  106  /ft 
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PRESSURE  RATIO,  p/p 


AXIAL  DISTANCE  FROM  NOZZLE  EXIT,  in.  (POSITIVE  DISTANCE  DOWNSTREAM) 


Fig.  29  Afterbody  Pressure  Distribution:  Comparison  Between  Experiment  and  Theory  for 
M„  =  1.2  and  Re  =  2.5  x  10® /ft 
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TABLE  I 

MEASUREMENT  UNCERTAINTY 


Parameter 

Uncertainty 

Boattail  Balance 

fn 

±0.  95  lbf 

fy 

±0.36  lbf 

fa 

±0.39  lbf 

Mm 

±1.  05  in. -lbf 

Afterbody  Balance 

% 

±7.  41  lbf 

fy 

±2.  79  lbf 

fa 

±0.24  lbf 

Mm 

±1.09  in. -lbf 

Boattail  Static  Pressures 

±1.  08  psf 

Nozzle  Exit  Static  Pressures 

±12.  96  psf 

Air  Venturi  Plenum  Pressure 

±1. 5  psi 

Ethylene  Venturi  Plenum 
Pressure 

±1. 5  psi 

Air  Venturi  Plenum  Temper¬ 
ature 

±1. 5°R 

Ethylene  Venturi  Plenum 
Temperature 

±1.  5°R 

Air  Venturi  Flow  Rate 

±0.  93  percent 

Ethylene  Venturi  Flow  Rate 

±0.  93  percent 
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TABLE  II 

NOZZLE  BO  ATT  AIL  SURFACE  PRESSURE  INSTRUMENTATION 


Model 
Station,  in. 

Angle  Orientation,  deg 

0 

45 

90 

135 

180 

270 

148. 101 

403 

412 

417 

426 

431 

440 

148. 811 

404 

418 

432 

441 

149.  591 

405 

413 

419 

427 

433 

442 

150.281 

406 

420 

434 

443 

150. 941 

407 

414 

421 

428 

435 

444 

151.  531 

408 

422 

436 

445 

152.071 

409 

415 

423 

429 

437 

446 

152. 531 

410 

424 

438 

447 

152. 951 

411 

416 

425 

430 

439 

448 

0  deg 


View  Looking  Upstream 
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TABLE  III 

FOREBODY  AND  AFTERBODY  SURFACE  PRESSURE  INSTRUMENTATION 


Model 
Station,  in. 

15.00 
29.  58 
41.  50 
54.  00 

j  66. 00 

78.  00 
92.  00 
103.  00 
110. 51 
114. 89 
119. 26 
123. 64 
128.01 
129.  20 

136. 97 
140. 22 
142. 72 
144.  72 
146.  22 
I  147.22 


180  deg 


View  Looking  Upstream 
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TABLE  IV 

BOUNDARY  CONDITIONS 


Region 

Boundary-  Condition 
Type 

Region  from  Which 
Boundary  Condition 

Is  Imposed 

A 

Solid  Boundary 

C 

Solid  Boundary 

E 

B 

Solid  Boundary 

D 

Solid  Boundary 

E 

C 

Pressure  Gradient 

- 

A 

D 

Pressure  Gradient 

B 

E 

Axial  Pressure  Gradient 

A  and  B 

Reference  Line 

Match  Pressure  (Transverse) 

between  A  and  B 

TABLE  V 

COMPUTATIONAL  METHODS 


Region 

Description 

Characteristics 

Method 

A 

Internal  Nozzle  Flow 
and  Plume 

Inviscid  -  Supersonic 

Lockheed  MOC 

B 

External  Flow 
(Isolated  Nozzle) 

Inviscid 

a.  Subsonic 

b.  Supersonic 

a.  Pratt  &  Whitney 
Transonic  Flow 
Program 

b.  Lockheed  MOC 

C 

Nozzle  Boundary 
Layer 

Viscous 

Patankar 

Option  I 

D 

External  Boundary 
Layer 

Viscous 

Patankar 

Option  I 

E 

Mixing  Layer 

Viscous 

Patankar 
Option  II 
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